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Abstract 

We consider a natural destruction process of an infinite recursive tree by removing each 
edge after an independent exponential time. The destruction up to time t is encoded by a 
partition 11(f) of N into blocks of connected vertices. Despite the lack of exchangeability, 
just like for an exchangeable fragmentation process, the process II is Markovian with 
transitions determined by a splitting rates measure r. However, somewhat surprisingly, r 
fails to fulfill the usual integrability condition for the dislocation measure of exchangeable 
fragmentations. We further observe that a time-dependent normalization enables us to 
define the weights of the blocks of n(t). We study the process of these weights and point 
at connections with Ornstein-Uhlenbeck type processes. 

Key words: Random recursive tree, destruction of graphs, fragmentation process, cluster 
sizes, Ornstein-Uhlenbeck type process. 


1 Introduction 

The purpose of this work is to investigate various aspects of a simple and natural fragmentation 
process on an infinite tree, which turns out to exhibit nonetheless some rather unexpected 
features. 

Specifically, we first construct a tree T with set of vertices N = {1,...} by incorporating 
vertices one after the other and uniformly at random. That is, 1 is the root, and for each 


*erich.baur@math.uzh.ch 

ljean.bertoin@math.uzh.ch 



vertex i > 2, we pick its parent rq according to the uniform distribution in {l,..,,i — 1}, 
independently of the other vertices. We call T an infinite (random) recursive tree. Recursive 
trees are especially useful in computer science where they arise as data structures; see e.g. the 
survey by Mahmoud and Smythe [14] for background. 

We next destroy T progressively by removing each edge e* connecting % to its parent Ui at time 
€i, where the sequence (e, : i > 2) consists of i.i.d. standard exponential variables, which are 
further independent of T. Panholzer [16] investigated costs related to this destruction process, 
whereas in a different direction, Goldschmidt and Martin [10] used it to provide a remarkable 
construction of the Bolthausen-Sznitman coalescent. We also refer to Kuba and Panholzer [13] 
for the study of a related algorithm for isolation of nodes, and to our survey [3] for further 
applications and many more references. 

Roughly speaking, we are interested here in the fragmentation process that results from the 
destruction. We represent the destruction of T up to time t by a partition IP(t) of N into 
blocks of connected vertices. In other words, if we view the fragmentation of T up to time t as 
a Bernoulli bond-percolation with parameter e - *, then the blocks of 11(f) are the percolation 
clusters. Clearly 11(f) gets finer as f increases, and it is easily seen from the fundamental 
splitting property of random recursive trees that the process II = (11(f) : f > 0) is Markovian. 
In this direction, we also recall that Aldous and Pitman [1] have considered a similar logging of 
the Continuum Random Tree that yields a notable fragmentation process, dual to the standard 
additive coalescent. We further point at the very recent work [11] in which the effects of repeated 
random removal of nodes (instead of edges) in a finite random recursive tree are analyzed. 

It turns out that II shares many features similar to homogeneous fragmentation processes 
as defined in [4, 5]. In particular, the transition kernels of II are very similar to those of a 
homogeneous fragmentation; they are entirely determined by the so-called splitting rates r, 
which define an infinite measure on the space of partitions of N. However, there are also 
major differences: exchangeability, which is a key requirement for homogeneous fragmentation 
processes, fails for n, and perhaps more notably, the splitting rates measure r does not fulfill 
the fundamental integral condition (6) which the splitting rates of homogeneous fragmentation 
processes have to satisfy. 

It is known from the work of Kingman [12] that exchangeability plays a fundamental role 
in the study of random partitions, and more precisely, it lies at the heart of the connection 
between exchangeable random partitions (which are discrete random variables), and random 
mass-partitions (which are continuous random variables). In particular, the distribution of an 
exchangeable random partition is determined by the law of the asymptotic frequencies of its 
blocks B , 

\B\ = lim n -1 #{i < n : i e B}. (1) 
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Even though n(£) is not exchangeable for t > 0, it is elementary to see that every block of 
II(t), say B(t), has an asymptotic frequency. However this asymptotic frequency is degenerate, 
\B(t)\ = 0 (note that if n were exchangeable, this would imply that all the blocks of n(t) would 
be singletons). We shall obtain a finer result and show that the limit 

lim n~ e *#{i < n : i G B{t)} (2) 

n—too 

exists in (0, oo) almost surely. We shall refer to the latter as the weight of the block B{t) 
(we stress that this definition depends on the time t at which the block is taken), and another 
natural question about the destruction of T is thus to describe the process X of the weights of 
the blocks of the partition-valued process n. 

Because n resembles homogeneous fragmentations, but with splitting rate measure r which 
does not fulfill the integral condition of the former, and because the notion (2) of the weight of 
a block depends on the time t, one might expect that X should be an example of a so-called 
compensated fragmentation which was recently introduced in [6]. Although this is not exactly 
the case, we shall see that X fulfills closely related properties. Using well-known connections 
between random recursive trees, Yule processes, and Polya urns, cf. [3], we shall derive a number 
of explicit results about its distribution. In particular, we shall show that upon a logarithmic 
transform, X can be viewed as a branching Ornstein-Uhlcnbeck process. 

The rest of this paper is organized as follows. In Section 2, we study the structure of the 
partition-valued process n which stems from the destruction of T, stressing the resemblances 
and the differences with exchangeable fragmentations. In Section 3, we observe that after a 
suitable renormalization that depends on t, the blocks of the partition n(t) possess a weight, 
and we relate the process of these weights to Ornstein-Uhlcnbeck type processes. 


2 Destruction of T and fragmentation of partitions 

The purpose of this section is to show that, despite the lack of exchangeability, the partition 
valued process n induced by the fragmentation of T can be analyzed much in the same way as 
a homogeneous fragmentation. We shall present the main features and merely sketch proofs, 
referring to Section 3.1 in [5] for details. 

We start by recalling that a partition n of N is a sequence ( 7 q : i G N) of pairwise disjoint 
blocks, indexed in the increasing order of their smallest elements, and such that Uj e N7q = N. We 
write V for the space of partitions of N, which is a compact hypermetric space when endowed 
with the distance 

d(7r,7r') = 1/ max{n e N : vr| [n] = 7rf [n] }, 
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where n \b denotes the restriction of tt to a subset B C N and [n] — {1,,.., n} is the set of the 
n first integers. The space Vb of partitions of B is defined similarly. 

We next introduce some spaces of functions on V. First, for every n > 1, we write D n for the 
space of functions / : V —» M which remain constant on balls with radius 1 /n, that is such that 
/(7r) = /(?/) whenever the restrictions 7T|[ n ] and r/|[ n ] of the partitions n and rj to [n] coincide. 
Plainly, D n C D n+1 , and we set 

Doo = 1J D n- 

n> 1 

Observe that D 00 is a dense subset of the space C(V) of continuous functions on V. 

In order to describe a family of transition kernels which appear naturally in this study, we 
first need some notation. For every block BCN, write B(j) for the j'-th smallest element of 
B (whenever it makes sense), and then, for every partition n G V, B o tt for the partition of B 
generated by the blocks B{rTi) = {B(j) : j e 7p} for i e N. In other words, B o tt is simply the 
partition of B induced by n when one enumerates the elements of B in their natural order. Of 
course, if the cardinality of B is finite, say equal to k e N, then B o tt does only depend on n 
through 7T|[fc], so that we may consider B o tt also for 7r e V[k\ (or tt e for any i > k). 

In the same vein, for partitions ij G Vb and every integer % > 1, we write rj on for the 

i 

partition of B that results from fragmenting the i-th block of rj by 7T, that is replacing the block 
rji in rj by rj t o tt. Again, if k — #B < oo, we may also take 7r e V\(\ for l>k. 

Finally, for every k > 2, we consider a random partition of N that arises from the following 
Polya urn. At the initial time, the urn contains k — 1 black balls labeled 1,..., k — 1 and a 
single red ball labeled k. Balls with labels k + 1, k + 2,... are colored black or red at random 
and then incorporated to the urn one after the other. More precisely, for n > k , the color given 
to the n + 1-th ball is that of a ball picked uniformly at random when the urn contains n balls. 
This yields a random binary partition of N; we write p*. for its law. We set 

OO 

r = Pa ’ ( 3 ) 
k =2 

which is thus an infinite measure on the set of binary partitions of N. 

Recall that each edge of T is deleted at an exponentially distributed random time, indepen¬ 
dently of the other edges. This induces, for every t > 0, a random partition 11(f) of N into 
blocks corresponding to the subsets of vertices which are still connected at time t. Observe 
that, by construction and the very definition of the distance on V, the process II has cacllag 
paths. 

We are now able to state the main result of this section. 
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Theorem 1 (i) The process II = (11(f) : t > 0) is Markovian and has the Feller property. We 
write G for its infinitesimal generator. 

(ii) For every n > 1, D n is invariant and therefore D^ is a core for G. 

(iii) For every f G D 00 and 77 E V, we have 

r (dvr)^ (f{r) ott) - . 

i 



We stress that this characterization of the law of the process fl is very close to that of a 
homogeneous fragmentation. Indeed, one can rephrase well-known results (cf. Section 3.1.2 
in [5]) on the latter as follows. Every homogeneous fragmentation process T = (T t : t > 0) 
is a Feller process on V, such that the sub-spaces D n are invariant (and hence D 0 0 is a core). 
Further, its infinitesimal generator A is given in the form 

= / s (d7r)^ (f(v°*) ~ /(h)) 

Jnev i 1 

for every / G D a 0 and 77 G V, where s is some exchangeable measure on V. More precisely, 
s({lpj}) = 0, where for every block B C N, 1b G Vb denotes the neutral partition which has a 
single non-empty block B, and 

s (| 7 T G V : 7 T|[ n ] 7 ^ l[n]}) < oo for all n > 2. (4) 

Observe that the measure r fails to be exchangeable, but it fulfills (4); indeed, one has 

n 

r ({ 7 T G V : 7 T| W ^ 1 [n] }) = 5^Pfc(P) =n- 1. 

k =2 



We shall now prepare the proof of Theorem 1. In this direction, it is convenient to introduce 
some further notation. Consider an arbitrary block BCN,a partition ij G Vb and a sequence 
7 rO = (7 : i G N) in V. We write 77 o V') for the partition of B whose family of blocks is given 
by those of 77 o for iGl In words, for each iGN, the i-th block of 77 is split according to 

i 

the partition 7 r^. Next, consider a probability measure q on V and a sequence ( 77 ^ : j g N) of 
i.i.d. random partitions with common law q. We associate to q a probability kernel Fr(-, q) on 
Vb, by denoting the distribution of 77 o by Fr (77, q) for every 77 G Vb- We point out that if 
q is exchangeable, then 77 o 7r^' ) has the same distribution as the random partition whose blocks 
are given by the restrictions tt^ of to 777 for i G N, and Fr(-,q) thus coincides with the 
fragmentation kernel that occurs for homogeneous fragmentations, see Definition 3.2 on page 
119 in [5]. Of course, the assumption of exchangeability is crucial for this identification to hold. 
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Note that the restriction of partitions of N to [n] is compatible with the fragmentation 
operator Fr(-, •), in the sense that 

(^° 7r( ' ) )|[n] =r 7lM 07r( ' } IN- (5) 

Proposition 1 For every n G N, the process II|[ n ] = (II|[ n ](t) : t > 0) obtained by restricting II 
to [n], is a continuous time Markov chain on V\ n ]. 

Its semigroup can be described as follows: for every s,t > 0, the conditional distribution of 
II|[ n ](s + 1) given II|[„](s) = g is Fr( 77 ,q t ), where denotes the distribution ofU(t). 

Proof: The case n — 1 is clear since IIj[i](t) = ({1}, 0, ...) for all times t > 0. Assume now 
n > 2. The proof relies crucially on the so-called splitting property of random recursive trees 
that we now recall (see, e.g. Section 2.2 of [3]). Given a subset BCN, the image of T by the 
map j i —> B(j) which enumerates the elements of B in the increasing order, is called a random 
recursive tree on B and denoted by T#. In particular, for B = [n], the restriction of T to the 
first n vertices is a random recursive tree on [n]. Imagine now that we remove k fixed edges 
(i.e. edges with given indices, say A,... ,4, where 2 < q < ... < i k < n) from T[ n j. Then, 
conditionally on the induced partition of [n], say g = (r/i,..., pk+i), the resulting k +1 subtrees 
are independent random recursive trees on their respective sets of vertices rjj, j — 1 ,..., k + 1. 

It follows easily from the lack of memory of the exponential distribution and the compatibility 
property (5) that the restricted process II|[ n ] = (ll|[ n ](*):* > 0 ) is a continuous time Markov 
chain. More precisely, the conditional distribution of II|[ n ](s + t ) given Il|[ n ](s) = n is Fr( 7 T, q t ), 
where Fr(-, q ( ) is here viewed as a probability kernel on Vwa. □ 

In order to describe the infinitesimal generator of the restricted processes ILr n ] for n e N, we 
consider its rates of jumps, which are dehned by 

r y = hm t _ 1 P(lI|[ n ](t) = 7 r), 
t^> o+ 

where now n denotes a generic partition of |n] which has at least two (non-empty) blocks. 
The rates of jumps r n determine the infinitesimal generator G n of the restricted chain II|[ n ], 
specihcally we have for / : Pr n i —>■ R and g G Pr n i 

Gnfiv) = (/( r ^° 7r ) - fivj) r n 

Tr£Vy n } i 

(recall that g o n denotes the partition that results from fragmenting the i-th block of g according 

i 

to 7 r). This determines the distribution of the restricted chain II|[ n ], and hence, letting n vary 
in N, also characterizes the law of II. Recall also that the measure r on V has been dehned by 
(3). 
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Proposition 2 For every n > 2 and every partition i r of [n] with at least two (non-empty) 
blocks, there is the identity 

r n = r(V n ), 

where = {rj G V : r/| [n ] = 7r}. 


Proof: This should be intuitively straightforward from the connection between the construc¬ 
tion of random recursive trees and the dynamics of Polya urns. Specifically, fix n > 2 and 
consider a partition n G P[ n ]. If n consists in three or more non-empty blocks, then we clearly 
have 

hmt- lp ( n| N (t) = tt) = 0, 

since at least two edges have to be removed from T|[ n ] in order to yield a partition with three or 
more blocks. Assume now that 7r is binary with non-empty blocks n\ and 7T2, and let k = min 7r 2 . 
Then only the removal of the edge may possibly induce the partition n, and more precisely, if 
we write r) for the random partition of [n] resulting from the removal of e*,, then the probability 
that r) = tt is precisely the probability that in a Polya urn containing initially k — 1 black balls 
labeled 1 ,..., k — 1 and a single red ball labeled k , after n — k steps, the red balls are exactly 
those with labels in 7 t 2 . Since the edge e*, is removed at unit rate, this gives 


lim t 1 P(n !M (/;) = 7r) = ry = p fc (P 7r ) 

i—>■ 04- 

in the notation of the statement. Note that the right-hand side can be also written as r(V n ), 
since peiV^) = 0 for all £ ^ k. □ 

Proposition 2 should be compared with Proposition 3.2 in [5]; we refer henceforth to r as 
the splitting rate of II. 

We have now all the ingredients necessary to establish Theorem 1. 

Proof of Theorem 1: From Proposition 1, we see that the transition semigroup of II is 
given by (Fr(-,q t ) : t > 0), and it is easily checked that the latter fulfills the Feller property; 
cf. Proposition 3.1 (i) in [5]. Point (ii) is immediate from the compatibility of restriction with 
the fragmentation operator, see (5). Concerning (iii), let / G D n . Since / is constant on 
{rj G V : r/|[ n ] = 7r}, it can naturally be restricted to a function / : V[ n ] R. By the 
compatibility property (5), with rv = r(TV) for a partition n' of [n] defined as in Proposition 2, 
we obtain 

/ r (dvr) {f(v ° “ f(vj) = (/falM ° 7 0 “ f(v\[n]j) = G n f(r}\ [n] ), 

Jn&v i 1 7v'GV [n] i 

where G n is the infinitesimal generator of the restricted chain II|[ n ] found above. This readily 
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yields (iii). □ 

Remark. It may be interesting to recall that the standard exponential law is invariant under 
the map t i-> — ln(l — e - *), and thus, if we set e* = — ln(l — exp(—e*)) (recall that e t is the 
instant at which the edge connecting the vertex % to its parent is removed), then (e,)j > 2 is a 
sequence of i.i.d. exponential variables. The time-reversal ft-)- — ln(l — e _t ) transforms the 
destruction process of T into a construction process of T defined as follows. At each time 
6j, we create an edge between i > 2 and its parent which is chosen uniformly at random in 
{l,...,i — 1}. It follows that the time-reversed process 11(f) = II(— ln(l — e~ 4 ) —), t > 0, 
is a binary coalescent process such that the rate at which two blocks, say B and B' with 
mini? < min B', merge, is given by ff{j G B : j < min B'\. This can be viewed as a duality 
relation between fragmentation and coalescent processes; see Dong, Goldschmidt and Martin [9] 
and references therein. 

Recall that for every k > 2, a random binary partition with law p^, resulting from the Polya 
urn construction possesses asymptotic frequencies in the sense of (1). This readily entails the 
following result. 


Proposition 3 For r-almost all binary partitions (Bi,B 2 ) G V, the blocks B\ and B 2 have 
asymptotic frequencies, and more precisely, we have 


f(\Bi\,\B 2 \)dr = / /( l-x,x)x 2 dx, 


i-p 


where f : [0, l] 2 —> M + denotes a generic measurable function. In particular, 


(1 — |£>i|)dr = 00 . 


>v 


Proof: Indeed, it is a well-known fact of Polya urns that for each k > 2, p/,.- almost every 
partition (Bi,B 2 ) has asymptotic frequencies with 1^1 + \B 2 \ = 1 and \B 2 \ has the beta 
distribution with parameters (1, k — 1), i.e. with density (k — 1)(1 — x) k ~ 2 on (0,1). Our claims 
follow immediately since 

OO 

y^(fc — 1)(1 — x) k ~ 2 = x~ 2 , x G (0,1). 

k =2 


□ 


It is interesting to recall that the splitting rates s of a homogeneous fragmentation must 
fulfill the integrability condition 



Ri|)ds < 00 , 
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which thus fails for r ! 

We next turn our attention to the Poissonian structure of the process n, which can be 
rephrased in terms similar to those in Section 3.1 of [5]. In this direction, we introduce a 
random point measure 

OO 

i=2 

on M + x V x N as follows. Recall that e t is the time at which the edge e t connecting the vertex 
% e N to its parent in T is removed. Immediately before time e*, the vertex i belongs to some 
block of the partition II(ej—), we denote the label of this block by 7j (recall that II is cacllag 
and that blocks of a partition are labeled in the increasing order of their smallest element). 
Removing the edge e* yields a partition of that block B = 11^ (e—) into two snb-blocks, which 
can be expressed (uniquely) in the form BoAj. This defines the binary partition A* and hence 
the point measure M unambiguously. The process II can be recovered from M, in a way similar 
to that explained on pages 97-98 in [5]. Roughly speaking, for every atom of M, say (t,A,k), 
II(£) results from partitioning the 7-tli block of II(t—) using A, that is by replacing Hk(t —) by 
Hk(t —) o A. Adapting the arguments of Section 3.1.3 in [5], we have the following result: 

Proposition 4 The random measure M is Poisson with intensity A <g) r <g) where A denotes 
the Lebesgue measure on R + and # the counting measure on N. 

Proof: Recall that we write lr n i = ([n], 0,...) for the partition of [n] which consists of a single 
non-empty block. Consider a Poisson random measure M' with intensity A ® r ® # as in the 
statement. Then M' has almost surely at most one atom in each fiber {t} < 8 > V <g) N, and the 
discussion below (4) shows that for each t' > 0 and every n e N, the number of atoms (t, n, k ) 
of M' with t < t r , 7 T|[ n ] jtz ij n j and k < n is finite. We may therefore define for fixed n G N 
a Pr n i-valued continuous time Markov chain (II'A (t) : t > 0) starting from II'^^O) = 1 [„.] as 
follows: If t is a time at which the fiber {t} xPxN carries an atom (t, n, k) of M' such that 
7 T|[ n ] 7 ^ 1[„] and k < n, then IT^t) results from Il'^^t—) by replacing its 7-th block n^(t—) 

by nj n] (t-) o7T| [n] . 

The sequence (IRkl (t) : n e N) is clearly compatible for every t > 0, in the sense that 
n ,[nI | M (t) = n'[ m l(t) = for all integers n > m. We deduce as in the proof of Lemma 3.3 in [5] 
that there exists a unique P-valued cacllag function (n'(t) : t > 0) such that II|r n ](t) = iT^t). 
Moreover, the i-th block II' (t) of II'(t) is given by the increasing union II' (t) = U ne pjII^(t), 
and it follows from the very construction of n't”! (t) that the process II' can be recovered from 
M' similarly to the description above the statement of the proposition. It remains to check that 
IT and II have the same law, which follows if we show that the restricted processes nW = II'^ 
and II |[„,] have the same law for each n G N. Fix n > 2, and denote by n a partition of [n] with 
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at least two non-empty blocks. From the Poissonian construction of IF^, with V n as in the 
statement of Proposition 2, we first see that 

lim r x p(n'M(t) = tt) = r(V n ). 

Next, if tt' 7 ^ it" G V[ n ], the jump rate of IF^ from tt' to tt" is non-zero only if tt" can be 
obtained from tt' by replacing one single block of tt' , say the k- th block Tr' kl by Tr' k o n, where 
7 T is some binary partition of [n]. This observation and the last display readily show that IF^I 
and n Hn] have the same generator, and hence their laws agree. □ 


3 The process of the weights 

Even though the splitting rates measure r of the fragmentation process II fails to fulfill the 
integral condition ( 6 ), we shall see that we can nonetheless define the weights of its blocks. The 
purpose of this section is to investigate the process of the weights as time passes. 


3.1 The weight of the first block as an O.U. type process 

In this section, we focus on the first block II 1 (t), that is the cluster at time t which contains 
the root 1 of T. The next statement gathers its key properties, and in particular stresses the 
connection with an Ornstein-Uhlcnbeck type process. 


Theorem 2 (i) For every t > 0, the following limit 

lim <n: j e Ili(i)} = Ad(t) 

n—too 

exists in (0, oo) a.s. The variable Ah(t) has the Mittag-Leffler distribution with parameter e~', 


P(AVt) G dx)/dx = - V ^- 4 - Tike-' + l^- 1 smlnke-')] 

TT k\ 

k =0 


equivalently, its Mellin transform is given by 

F(g + 1) 


E (V?«) = 


r(e-‘9+l)’ 


q > 0. 


(ii) The process {X\ (t) : t > 0) is Markovian; its semigroup Pt(x, •) is given by 

P t (x,-)=F(x e ~ t X 1 (t)e-). 
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(iii) The process Y(t ) = lnAd(t), t > 0, is of Ornstein-Uhlenbeck type. More precisely, 


L(t ) = Y(t) + / y(s)ds, t > 0, 

Jo 

is a spectrally negative Levy process with cumulant-generating function 


K,(q) — lnE(exp(gL(l))), q > 0, 


given by 


«(e) = #(<? + 1 ), 


where ip denotes the digamma function, that is the logarithmic derivative of the gamma function. 


In the sequel, we shall refer to A \{t) as the weight of the first block (or the root-cluster) at 
time t. Before tackling the proof of Theorem 2, we make a couple of comments. 

Firstly, observe from (i) that linp^oo E(A^ 1 (t) <? ) = T(q+ 1), so that as t —> oo, Y(t) converges 
in distribution to the logarithm of a standard exponential variable. On the other hand, it is well- 
known that the weak limit at oo of an Ornstein-Uhlenbeck type process is self-decomposable; 
cf. Section 17 in Sato [18]. So (iii) enables us to recover the fact that the log-exponential 
distribution is self-decomposable; see Shanbhag and Sreehari [19]. 

Secondly, note that the Levy-Khintchin formula for k reads 

r0 

«(?) = -7? + / (e gx - 1 - qx) T --ydx, 

J- oo (! - e ) 

where 7 = 0.57721... is the Euler-Mascheroni constant. Indeed, this follows readily from the 
classical identity for the digamma function 

r 1 1 — x i 

ip(q + 1) = -7 + / -- dx. 

Jo 1 - x 

In turn, this enables us to identify the Levy measure of L as 

A(dx) = e x (l — e x ) _ 2 da;, x E (— 00, 0 ). 

Since the jumps of L and of Y coincide, the Levy-Ito decomposition entails that the jump process 
of Y = In Ad is a Poisson point process with characteristic measure A. In this direction, recall 
from Proposition 3 that the distribution of the asymptotic frequency of the first block under 
the measure r of the splitting rates of II is (1 — y)~ 2 dy, y G (0,1), and observe that the image 
of the latter by the map y 1 —y In y is precisely A. This should of course not come as a surprise. 
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We shall present two proofs of Theorem 2(i); the first relies on the well-known connection 
between random recursive trees and Yule processes and is based on arguments due to Pitman. 
Indeed, < n : j G Ili(f)} can be interpreted in terms of the two-type population system 
considered in Section 3.4 of [17], as the number of novel individuals at time t when the birth 
rate of novel offspring per novel individual is given by a = e _t , and conditioned that there are 
n individuals in total in the population system at time t. Part (i) of the theorem then readily 
follows from Proposition 3.14 in connection with Corollary 3.15 and Theorem 3.8 in [17]. For 
the reader’s convenience, let us nonetheless give a self-contained proof which is specialized to 
onr situation. We further stress that variations of this argument will be used in the proofs of 
Proposition 5 and Corollary 2. 

First proof of Theorem 2(i): Consider a population model started from a single ancestor, 
in which each individual gives birth a new child at rate one (in continuous time). If the ancestor 
receives the label 1 and the next individuals are labeled 2, 3,... according to the order of their 
birth times, then the genealogical tree of the entire population is a version of T. Further, 
if we write Z(s) for the number of individuals in the population at time s, then the process 
( Z(s ) : s > 0) is a Ynle process, that is a pure birth process with birth rate n when the 
population has size n. Moreover, it is readily seen that the Yule process Z and the genealogical 
tree T are independent. 

It is well-known that 

lim e~ s Z(s) = W almost surely, 

s—>• oo 

where W has the standard exponential distribution. As a consequence, if we write r n = inf{s > 
0 : Z(s) — n} for the birth-time of the individual with label n, then 

lim ne~ Tn = W almost surely. (7) 

n—>• oo 

Now we incorporate destruction of edges to this population model by killing merciless each 
new-born child with probability 1—p E (0,1), independently of the other children. The resulting 
population model is again a Yule process, say Z^ = (Z^(s) : s > 0), but now the rate of 
birth per individual is p. Therefore, we have also 

lim e~ ps Z ( ' p \s) = W (p) almost surely, 

s—>• OO 

where W^ is another standard exponential variable. We stress that is of course correlated 
to W and not independent of T, in contrast to W. 

In this framework, we identify for p = e~ t 

#{j<n:jen 1 (i)} = %) 
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and therefore 


lim e pTn #{j < n : j G IR^)} = almost surely. 

n—>-oo 

Combining with (7), we arrive at 

lim n p #{j < n : j G II 1 (f)} = almost surely, 

which proves the hrst part of (i). 

Now recall that the left-hand side above only depends on the genealogical tree T and the 
exponential random variables e t attached to its edges. Therefore, it is independent of the Yule 
process Z and a fortiori of W. Since both W and are standard exponentials, the second 
part of (i) now follows from the moments of exponential random variables. □ 

The second proof of Theorem 2(i) relies on more advanced features on the destruction of 
random recursive trees and Poisson-Dirichlet partitions. 

Second proof of Theorem 2(i): It is known from the work of Goldschmidt and Martin [10] 
that the destruction of T bears deep connections to the Bolthausen-Sznitman coalescent. In 
this setting, the quantity 

#{j < n: j e ni(t)} 

can be viewed as the number of blocks at time t in a Bolthausen-Sznitman coalescent on 
[n] — {1,..., n} started from the partition into singletons. On the other hand, it is known that 
the latter is a so-called (e - *, 0) partition; see Section 3.2 and Theorem 5.19 in Pitman [17]. Our 
claims now follow from Theorem 3.8 in [17]. □ 

Proof of Theorem 2 (ii): Let IT, be an independent copy of the process IR. Fix s,t > 0 and 
put B — II 1 (s), C — IT, (t). Recall that B(j ) denotes the j-th smallest element of B, and B{C) 
stands for the block { B(j ) : j e C}. By Proposition 1, there is the equality in distribution 
II 1 (s + 1) — B(C). From (i) we deduce that 

B(n) ~ (n/Xi(s)) e almost surely as n oo, 

and similarly C(n) ~ ( n/X[(t)) e as n —> oo, where X[(t) has the same law as Xi(t) and is 
further independent of (Ad(r) : r > 0). It follows that there are the identities 
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A"i(s + t) — lim m e ( + °#{j < m : j E Ui(s + t)} 

m—> oo 

= lim ((B(C)(n))~ e (s+t) n') 

n—>oo V / 

= lim ("(fi(C(n))) _e (s+t) n') 

n—>-oo \ / 

= xr t ( S )A-(t). 


Here, in the next to last equality we have used the fact that the n-th smallest element of B(C) 
is given by the C(n)- th smallest element of B, and for the last equality we have plugged in the 
asymptotic expressions for B{n) and C(n) that we found above. Our claim now follows easily. 
□ 

We point out that, alternatively, the Markov property of X l can also be derived from the 
interpretation of #{j < n : j E n^f)} as the number of blocks at time t in a Bolthausen- 
Sznitman coalescent on [n]; see the second proof of Theorem 2(i) above. 

Proof of Theorem 2(iii) : We first observe from (ii) that the process Y is Markovian with 
semigroup Qt(y,-) given by 

Q t (y,') = ne~ t y + Y(t)e-). 

Next, recall from the last remark made after Theorem 2 that the function q i —> n(q) = q'4>{q + 1) 
is the cumulant-generating function of a spectrally negative Levy process, say L = ( L(t) : t > 0). 
Consider then the Ornstein-Uhlenbeck type process U = {U(t) : t > 0) that solves the stochastic 
differential equation 

U(t) = L{t) - [ U(s)ds, 

Jo 

that is, equivalently, U(t) = e - * e s dL(s). Then U is also Markovian with semigroup Rt(u, •) 
given by 

Rt(u, ■) = P(e _t « + U(t) e ■). 

So to check that the processes Y and U have the same law, it suffices to verify that they have 
the same one-dimensional distribution. 

The calculations of Section 17 in Sato [18] (see Equation (17.4) and Lemma 17.1 there) show 
that for every q > 0, 

E {exp{qU(t))) = exp (J hi(e~ s q)ds 

Now observe that 



Kle 


S q)ds 



e s g^(e s q + l)ds = lnT(g + 1) — lnT(e t q+ 1), 
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and hence 


E (exp(„t/(t))) = ^, + 1> = E (exp(«y(())), 

T(e tq + 1) 

where the second identity stems from Theorem 2(i). □ 

Remark. It would be interesting to understand whether with probability 1, the block II 1 (t) 
has weights in the sense of Theorem 2(i) simultaneously for all t > 0 (note that the asymptotic 
frequencies are equal to zero on (0, oo) a.s.), and whether t H» X\ (t) is cacllag. In this direction, 
we recall that each block of a standard homogeneous fragmentation process possesses asymptotic 
frequencies simultaneously for all t > 0 a.s., see Proposition 3.6 in [5]. Moreover, if B denotes a 
block of such a process, then the process t (->• \B(t)\ is cadlag. Here it should be observed that 
the first block Hi is the only block which is decreasing, in the sense that Hi(f') is contained 
in Ih(t) a.s. whenever t! > t. This however does not imply monotonicity of X\(t) in t, which 
is crucial ingredient for the proof of the mentioned properties in the case of a homogeneous 
fragmentation. 

3.2 Fragmentation of weights as a branching O.U. process 

We next turn our interest to the other blocks of the partition n(f); we shall see that they also 
have a weight, in the same sense as for the first block. In this direction, it is convenient to 
write first Tj for the subtree of T rooted at i > 1: in particular Ti = T. Then for t > 0, we 
write Ti(t) the subtree of Tj consisting of vertices j G Tj which are still connected to i after 
the edges e*, with ej~ < t have been removed. Note that for i > 2, Tj(f) is a cluster at time t if 
and only if 6j < t, an event which has always probability 1 — e~ f and is further independent of 
7j(f). On that event, the vertex set of Ti(t) is a block of the partition n(f), and all the blocks 
of n(f) arise in this form. 

Lemma 1 For every t > 0 and i € N, the following limit 

lim < n : j £ T.ft)} = p^t ) 

n—>• oo 

exists in (0, oo) a.s. Moreover, the process 


Pi = (. pi(t ) : t > 0) 


has the same law as 

(PfX, (t) : t > 0), 


where fa denotes a beta variable with parameter (1 ,i — 1) and is further independent of X x (t). 
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In particular, the positive moments of pi(t) are given by 


E(P?W) 


r(g + i)r(i) 

T(e~ t q + i) 


q > 0. 


Proof: The recursive construction of T and Tj has the same dynamics as a Polya urn, and basic 
properties of the latter entail that the proportion /3j of vertices in Tj has the beta distribution 
with parameter (1, % — 1). Further, enumerating the vertices of Tj turns the latter into a random 
recursive tree. Our claim then follows readily from Theorem 2. □ 

Lemma 1 entails that for every j 6 N, the i-th block flj(f) of 11(f) has a weight in the sense 
of (2), a.s. We write Xfit) for the latter and set X(t) = (X\{t), X- 2 (t),.. .). We now investigate 
the process X = (X(t) : t > 0). 

Firstly, using the functional equation of the gamma function, the integral representation of 
the beta function and the expression for the moments of X b(t) from Theorem 2(i), an easy 
calculation shows 


E 



E (X q (t)) + E 


oo 

X 

, i=2 




to-i) r(g) 
(e _t g — 1) r(e _t g) ’ 


( 8 ) 


provided q > e*. In particular, for every t > 0, the Xfit) can be sorted in the decreasing order. 
We write Xl(f) for the sequence obtained from X(t) by ranking the weights X { (t) decreasingly, 
where as usual elements are repeated according to their multiplicity. For q > 0, let 


i q ^ — < x = (xi,...) : Xi > x 2 > ■ ■ ■ > 0, and ^://( < 


oo 


1=1 


endowed with the ^-distance. Similarly, denote by the space of ordered sequences of 
positive reals, endowed with the f' 00 -distance. For the process X^, we obtain the following 
characterization. 


Corollary 1 Let T G (0, oo], and set q = e T (with the convention e°° = oo). Then the process 
X^ = (X'*'(t) : t < T) takes its values in i q ^ and is Markovian. More specifically, its semigroup 
can be described as follows. For s,t > 0 with s + t < T, the law of X^s + t) conditioned on 
Xf(s) = (an,...) is given by the distribution of the decreasing rearrangement of the sequence 
(x® t x^ > : i,j G N), where ((x^ ,...) : % G N) is a sequence of independent random elements in 
£ q each of them distributed as Xl(i). 

Proof: The fact that X^(f) G £ q ' for t < T follows from (8). The specific form of the semigroup 
can be deduced from Theorem 1 and a variation of the arguments leading to Theorem 2(ii). □ 
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We now draw our attention to the logarithms of the weights. In this direction, recall that 
for a homogeneous fragmentation, the process of the asymptotic frequencies of the blocks bears 
close connections with branching random walks. More precisely, the random point process with 
atoms at the logarithm of the asymptotic frequencies and observed, say at integer times, is 
a branching random walk; see [8] and references therein. This means that at each step, each 
atom, say y , is replaced by a random cloud of atoms located at y + z for z G Z, independently 
of the other atoms, and where the random point process Z has a fixed distribution which 
does not depend on y nor on the step. In the same vein, we also point out that recently, 
a natural extension of homogeneous fragmentations, called compensated fragmentations, has 
been constructed in [7], and bears a similar connection with branching Levy processes. Note 
that compensated fragmentations are valued processes, in contrast to XL 

Our observations incite us to introduce 

y t = 

Theorem 3 The process with values 
Ornstein- Uhlenbeck branching process, 
abilities can be described as follows: 

For every s,t > 0, the cojiditional law 

of 

where the point measures 

3 =i 

are independent and each has the same law as 34- 

Furthermore, the mean intensity of y t is determined by 


Xi(t), t> 0. 


i =1 


in the space of point measure y = (y t : t > 0) is an 
in the sense that it is Markovian and its transition prob- 

of y s+t given 34 = ^y,. * s given by the distribution 

CO CO 

^e-*j/i+cO 

i= 1 j =1 


CO 




Proof: The hrst claim follows from Theorem 1 and an adaption of the argument for Theo¬ 
rem 2(iii). For the second, we note that 


E 


e qy y t (dy) 
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and the right hand side has already been evaluated in (8). □ 

Next we give a description of the finite dimensional laws of X(f) = (Ah(f), X 2 (t), ...) for 
t > 0 fixed. In this direction, it is convenient to define two families of probability distributions. 

The first family is indexed by j G N and t > 0, and is defined as 

nAk)=( k k _~ 2 _\e-r- i - 1 (l-e- t y i k>j + 1. 

Note that the shifted distribution p,j t t(k) = Hj,t(k + 1), k > j, is sometimes called the negative 
binomial distribution with parameters j and 1—e _t , that is the law of the number of independent 
trials for j successes when the success probability is given by 1 — e“b 

The second family is indexed by j G N and k > j, and can be described as follows. We 
denote by 9j k the probability measure on the discrete simplex A k j = {(k\,..., kj) G kb : 
ki H— • + kj — k}, such that 9j^ (k \,..., kj) is the probability that on a random recursive tree 
of size k (that is on a random tree distributed as T| [*,]), after j — 1 of its edges chosen uniformly 
at random have been removed, the sequence of the sizes of the j subtrees, ordered according to 
the label of their root vertex, is given by (ki ,..., kj). 

Remark. The distribution 9j^ is equal to 5k for j — 1. For j = 2, Meir and Moon [15] found 
the expression 

k 

9 2 ,k(ki , k 2 ) = -r —.——ttt--T, h, k 2 e N with k l + k 2 = k, 

k 2 {k 2 + 1 ){k - 1) 

with 9 2i k{k\,k 2 ) = 0 for all other pairs (ki,k 2 ). Generalizing the proof of this formula given 
in [15] to higher j, we find (k > j > 3 and k\ + • • • + kj = k) 


9j,k(ki, k 2 ,..., kj) 


(ki — 1)!(^2 — 1 )! • • • (kj — 1 )! fk — ij\ 

(k -l)\(k -l)---(k-(j-l)) kj ) 


(k—kj—kj-i)/\(£j—l) 

- E 

£,- 1 = 3-2 


k — kj — lj-\ 
kj -i 


(fc-Ef=2 3-1) 

x •" x E 
< 2=1 


k - EL 3 k > 

k 2 


Proposition 5 Let j G N, q ly . .., 1 > 0, and set kj + i = 1. The Mellin transform of the 

vector (Ah(f),... A j + i(t)) for fixed t > 0 is given by 


E (A'?' (()■■. X]Z‘ (()) = V Ht (k) V (*,.fcj) 


m 


3+1 

n 


k=j +1 


fci H-ffe 7 -=fc—1 


r(® + h 

T(ge _t + /c) 11 T(fci) 


Remark. By plugging in the dehnition of Hjj(k), one checks that the right hand side is finite. 
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Proof: Fix t > 0, and set p = e~*. For ease of notation, we write 11* and X* instead of H;(t) 
and Xi(t). Furthermore, fix an integer j G N and numbers ki,... ,kj € N. For convenience, set 
k — ki + ■■■ + kj + kj + 1 , with kj+ 1 = 1. We first work conditionally on the event 

Ak u ...,kj = { minll i+ i = k, #(fli n [k]) = h,.. ., #(II i n [k]) = kj). 


We shall adapt the first proof of Theorem 2(i). Here, we consider a multi-type Yule process 
starting from k individuals in total such that k, of them are of type i, for each % — 1 ,,,., j + 1 . 
The individuals reproduce independently of each other at unit rate, and each child individual 
adopts the type of its parent. Then, if Z(s) stands for the total number of individuals at time 
s, we have that lim^oo e~ s Z(s) = 7 (k) almost surely, where y(/c) is distributed as the sum of k 
standard exponentials, i.e. follows the gamma law with parameters (k, 1). Now assume again 
that each new-born child is killed with probability l-p 6 ( 0 , 1 ), independently of each other. 
Writing Z^ ,p \s) for the size of the population of type i at time s (with killing), we obtain 


lim e~ ps Z^ p \s) = 7 (i - p) (fci), i = l,...,j + l, 

S —^OO 


where the 7 ^ ,p \ki) are independent gamma( ki, 1 ) random variables (they are however clearly 
correlated to the asymptotic total population size 7 (k)). From the arguments given in the first 
proof of Theorem 2(i) it should be plain that conditionally on the event ^ we have for 
the weights X t the representation 


7 ^ p \ki) 
7 p {k) 


i = l,---,j + 1 , 


and the X* are independent of 7 (k). Now let 51 ,..., q 3 +\ > 0 and put q = q\ H— • + q 3 +\- Using 
the expression for the X* and independence, we calculate 


E (7(*n E (AT ■ • • A’7 j A h . tj ) = n r( ‘Yf‘ ) 

i= 1 ^ l ' 


(9) 


Therefore, again with k ]+ \ = 1, 

OO 

E(AT---.yy)= v 


k=j +1 


r (k) -p-r r(g*; + ki 

bYW fj r(fci) 

fclH- \-kj=k—l 




With k — k\ + ■■■ + kj + 1 as above, we express the probability of A^.^kj as 


P (Ak!,...,^) = p (#(n x n [k]) = ki,..., #(n i n [k]) = kj \ minn j+ i = k) P (minn i+ i = k ). 
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By induction on j we easily deduce that min n,- + i — 1 is distributed as the sum of j independent 
geometric random variables with success probability 1 —p, i.e. minn ?+1 — 1 counts the number 
of trials for j successes, so that P (ruin ITj+i = k) = )• Moreover, it follows from the very 

definition of the blocks II* and the fact that the exponentials attached to the edges of TVi are 
i.i.d. that 


p (#(iii n [k]) = ki,. .., #(n i n [k]) = % | min n j+1 = k) = e jtk -i{k u ..., %). 

This proves the proposition. □ 

We finally look closer at the joint moments of X\ (t) and X 2 (t) when t tends to zero. We 
observe 9 1<k = S k and = (e _ *) fc_2 (l — e _t ), so that 

E(AT(i)A'? W ) = (1 -e-*)r(® + 1 )p(k - 

Now assume > 1- From the last display we get 

Urn i E (AT(t)A? (*)) = f; Tfa+orfe + i) 
t-+0+t “ T(qi + q 2 + k + l) 


OO »1 

^ k / (1 — x) <?1+fc_1 a;' ?2 da; = / (1 — a;) 91 aP 2-2 da;, 

fc=i ^ 


which one could have already guessed from Proposition 3. 


3.3 Asymptotic behaviors 

We shall finally present some asymptotic properties of the process X of the weights. To start 
with, we consider the large time behavior. 

Corollary 2 As t —> oo, there is the weak convergence 

(Xi(t) : i e N) =► (Wi : i G N), 

where on the right-hand side, the W.\ are i.i.d. standard exponential variables. 

Remark. This result is a little bit surprising, as obviously II(oo) is the partition into singletons. 
That is Ilj(oo) is reduced to {i} and hence has weight 1 if we apply (2) for t — oo. In other 
words, the limits n —> oo and t —> oo may not be interchanged. 
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Proof: Fix j e N arbitrarily and consider the event A(t) = {minllj + i(£) — j + 1}. Recall 
from the proof of Proposition 5 that P(A(t)) = / M,j(j +1) and note that this quantity converges 
to 1 as t —» oo. Further, on the event Aft), we have also bh(t)n[j + l] = {i} for all 1 < i < j + 1, 
that is Aft) = again in the notation of the proof of Proposition 5. 

Take q., ..., q j+1 > 0. Applying (9), we get 


j+i 

| AW) =n r («f + l). 

%— 1 

which entails that (Xi(t),..., X ]+ \(t)) converge in distribution as t —» oo towards a sequence 
of j + 1 independent exponentially distributed random variables. □ 

We next consider for t > 0 fixed the behavior of X n (t) as n — * oo. 

Corollary 3 Let t > 0. As n —* oo, there is the weak convergence 

rf* X n (t) =► P e_t Xi (t), 

where V denotes an exponential random variable of parameter (1 — e - *) -1 which is independent 
ofX.it). 

Proof: In the notation of Lemma 1, we have X n (t) = Pi( ni t){t ), where (again with the conven¬ 
tion ci = 0) 

i(n, t ) = min < j > 1 : ^ = n 

l i= i 

From Lemma 1 we know that X n (t) = fipn ^X\{t) in distribution, where i(n,t) and Pi{n,t) are 
both independent of X\ (t). By the law of large numbers, i(n,t) ~ n( 1 — e - *)^ 1 almost surely. 
Writing 

A(n,t) = i(n,i) _1 

and using the fact that kf3k converges in distribution to a standard exponential random variable 
as k —y oo, the claim follows. □ 

We finally look at the case t —> 0+. From Theorem 2(i) it follows that 

lim Xfit) = 1 in probability. (10) 

In fact, Xi(t) tends to be the largest element of the sequence X when t approaches zero. 

Corollary 4 

lim P (W(t) > X t (t) for all i > 2) = 1. 
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Proof: We consider the complementary event. We have 


P (Xi(t) > Xi(t) for some i > 2) < P (A/(f) > 1/2 for some i > 2) + P (Xi(t) < 1/2), 

and the second probability on the right hand side converges to zero as t —> 0+ by (10). For the 
first probability, we have by Lemma 1, with /3j denoting a beta(l,i — 1) random variable, 

°o 

F(Xi(t) > 1/2 for some i> 2) < (1 - e _t ) J^P [piXf (t) > (l/2) e ^ . 

i =2 

Using independence and the expression for the moments of X\(t) from Theorem 1, we obtain 
for t > 0 such that e* < 2, E(/3iXf (t)) < 2/i and Vai(/3iXf (f)) < 24 /i 2 . Therefore, for such t 
and i > 10, by Chebycheff’s inequality, with C = 10 4 , 

P (iftjff’t i) - E(AA'f(i))| > (1/2)'' - 2/,:) < 2 

This shows 

P {Xi(t) > 1/2 for some i > 2) < (1 — e -< ) 

□ 

Concluding remark. As it should be plain from the introduction, the set of vertices 

cf’m = 0 <11 : j 6 n,(«)} 

form the percolation clusters of a Bernoulli bond percolation with parameter p = e - * on a 
random recursive tree of size n on the vertex set {l,...,n}. Our results of Section 3 can 
therefore be understood as results on the asymptotic sizes of these clusters when n tends to 
infinity. 

Cluster sizes of random recursive trees were already studied in [6] when the percolation 
parameter p satisfies p = p(n) = 1 —s/ In n+o(l/ In n) for s > 0 fixed. The analysis was extended 
in [2] to all regimes p(n) —> 1. It shows that in these regimes, the root cluster containing 1 has 
always the asymptotic size ~ n p ( n \ while the next largest cluster sizes, normalized by a factor 
(1 — p(n))~ 1 n~ p ^ n \ are in the limit given by the (ranked) atoms of a Poisson random measure 
on (0, oo) with intensity a -2 da. 

The regime of constant parameter p = e - * considered here deserves the name “critical”, since 
in this regime, the root cluster and the next largest clusters have the same order of magnitude, 
namely n p . This is already apparent from Lemma 1. More precisely, Corollary 2 readily shows 


uxj 1 

+ cEhU°( ( ) 
*=10 1 


as t —^ 0 T . 
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that in fact 


limliminf P (there exists i > 2 such that #C,- n) (f) > #C'{ n ' ) (f)j = 1. 

tfoo n—too \ J 
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